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We calculate a grand partition function of the attractive Bose gas in the 
j~^ ' infinite space within some approximations. Using the idea of the Yang-Lee 

On ' zeros, it is proved that the gas-liquid condensation occurs before the conven- 
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tional condition of the Bose-Einstein condensation is satisfied. Further, it is 



C^ ' pointed out that Bosons with a zero momentum play a role of a trigger to 

a" 

this gas-liquid condensation. We discuss its implication to the trapped atomic 



gas. 
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The study of the relationship between the Bose-Einstein condensation (BEC) and the 
gas-hquid condensation (GLC) is a long-standing problem |I|. Recently the experimental 
realization of BEC in dilute atomic vapors changed this academic problem to a realistic one 
0. Normally, the BEC, a condensation into a lowest energy level at low temperature and 
high density, is thought to be an essentially different phenomenon from the GLC in following 
points: (1) The BEC is caused by the Bose statistics, not by the interparticle interaction as 
in the GLC. (2) The BEC is sometimes described as a condensation in momentum space, 
while the GLC occurs in coordinate space. 

An interesting point is that the reason by which we distinguish them is not so obvious 
as it looks. Both condensations have a similar thermodynamic manifestation as a first-order 
phase transition. Further, in two types of the quantum gas, the attractive Fermi gas and the 
attractive Bose gas, the relationship between the BEC and the GLC manifests itself quite 
differently. 

For Fermions, the circumstance with low temperature and high density stabilizes the 
Cooper pairs: a BEC in a general sense [Q. Because of the Fermi statistics, however, 
two Fermions experience a strong repulsive force in the short distance. Hence, when the 
attractive force is increased within the Bardeen-Cooper-Schriefer (BCS) model, the GLC is 
impossible, a rigorous proof of which is given recently @]. 

For Bosons, they remain in the gas state at high temperature and low density, because 
the quantum statistics plays a minor role. At low temperature, however, the Bosons have no 
large positive kinetic energy to stabilize the system behavior, especially in the BEC state. 
Thus, the attractive force plays a dominant role, so that the compressibility is no longer 
positive definite. The dilute Bose system will collapse into the dense one, leading to the 
GLC. Conversely, when the particle density increases by the GLC, an overlapping of the 
wave function is likely to cause the BEC. 

There are some kinematical or thermodjTiamical evidences suggesting the instability of 
the attractive Bose gas. 

(1) The Bogoliubov model says that if the interaction between the Bosons is attractive. 



the velocity of sound propagating on the Bose-Einstein condensate would be imaginary, 
corresponding to a divergence of the density fluctuation and leading to a drastic change of 
the whole system p. 

(2)A variational argument of the many-Boson system shows that a many-body self- 
binding state (liquid) is more likely to be stable than a gas of the bound dimers [Q . 

(3) Recent stability analysis of the Gross-Pitaevskii equation in the trapped atomic gas 
with the attractive interaction indicates a collapse of the BEG at the thermal equilibrium. 

These features suggest that an another type of relationship between the BEG and the 
GLG exists in the attractive Bose gas. The GLG in the classical imperfect gas is a long- 
standing and difficult problem [0. In the attractive Bose gas, however, there is a clear 
physical reason for the instability originated from the Bose statistics. In view of this, at 
least for the instability mechanism to the liquid, a model of the GLG in the attractive Bose 
gas would be simpler than that in the imperfect classical gas or in the attractive Fermi gas. 
With this in mind, it seems quite natural to consider the BEG and GLG on a common 
ground. 

In this paper, we consider a spinless Bose gas with a repulsive core represented by H^e 
and a weak attractive s-wave pairing interaction g{< 0) |^: 

H = Y^ eka\ak + Hre + ^Yl ala-kO'-k'ak'. (1) 

it ^ k,k' 

(In Eq.(l), Hre is included only to assure universal natures of the many-body system, so 
that this paper does not deal with it explicitly.) We assume that the Hre and the attrac- 
tive interaction makes an additive contribution to the free energy because of their differ- 
ent effective ranges, with a result that the grand partition function Z{fi) is factorized as 
Zo{fi)Zre{fi')Zat{fi')- Further we assume that diluteness of the system allows a contact inter- 
action in f7 as a first approximation. What concerns us is a singularity induced by Zat{lj)- 
We point out the following: (i) When the temperature decreases and the density increases, 
the grand partition function of the attractive Bose gas becomes zero at a negative critical 
value of the chemical potential ^c{< 0) . (ii) The Bosons with a zero momentum play a 



special role in this instability. 

The GLC can be considered as a singularity in the isothermal p — V diagram. The 
pressure p and the density p are given by, 

7^ = li- ^> (2) 

^ = hm A f l!^") , (3) 

where Zy is the grand partition function in the volume V and /i is a chemical potential. A 
zero of ^(yu) defines the GLC, giving the isothermal line a continuous but not differentiable 
point (Yang-Lee's zero 0). On the contrary, a pole of 2'(/i) defines the BEG. Because of 
the logarithmic function in Eqs.(2) and (3), the pole gives the isotherm near the BEG point 
a similar shape to that by the zero. In contrast to the GLG, however, the isotherm is not 
only continuous but also differentiable at the BEG point of the free Bose gas. 

The 2^(/i) which explains the GLG should have a following feature. Generally in the 
condensed matter, local changes induced by an external perturbation, or local interactions 
responsible for the thermodynamic quantities in the normal phase are well described by a 
sum of the disconnected ring diagram like Fig. 1(a). To describe a macroscopic change like the 
GLG, however, it is essential to include much greater networks of the interaction extending 
to all particles which consist of the system. As the particle number increases, a variety 
of such globally connected diagrams increases rapidly. Further, this connected networks is 
necessary for satisfying the quantum statistics correctly. When two particles having a same 
momentum and belonging to different ring diagrams are exchanged, a connected graph is 
produced as a result (for Fig. 1(a), 1(b) is produced). Hence, to include a symmetrized or 
antisymmetrized state into calculation, the sum of the disconnected ring diagrams must be 
supplemented by the corresponding connected diagram. Summing up all these diagrams for 
Z(/x) is an essential step to describe the macroscopic change caused by the quantum statistics 
such as the GLG in the quantum system. (This viewpoint was originally formulated in the 
case of the attractive Fermi gas: superconductivity [|lOl |1TT| .) 



For completeness, we repeat the formalism in Ref. [|T0| |]TT| in a simplified form, and apply 
it to the Boson system. The connected diagram like Fig. 1(b) is made of a combination of 
elementary diagrams like Fig. 1(c) consisting of 2s particle lines. Since each graph in Fig. 1(c) 
consists of Boson lines with a single {l,p), a sum over {l,p), 



is an elementary unit |]T2[. (Since each interaction line is connected to two Ks, each Ks has 
{g/V)^.) Only after summing the -ftTs's over all ways of connecting them by the interaction 
lines, we get an expected Z{fi) . 

Consider a connected diagram like Fig. 1(b) where the Ks appears z/^ times as in {z/s} = 
(z/i, z/2, . . .) (In Fig.l(b), {z/^} = (3, 3, 0, . . .)), and the Kg^s are connected to each other by n 
interaction lines. 

(1) There are Ugl ways of rearrangement which leaves the diagram invariant. (2) Since 
there are a number of ways of distributing frequency / and momentum p to each Kg, each 
interaction line connecting the -ft's's has an individuality characterized by {l,p) and {l',p') 
of particles which enter or emerge at both ends. This allows us n\ ways of rearrangements 
which produce different diagrams. 

With this in mind, we get, 

where 2s in the denominator is a number of rotations which leave the Kg invariant. 
If the sum over Ug can be carried out independently to n, we simply obtain — ^-^ — 



ZoZr 



nlll^ exp{—^), but in reality they are related to each other by n = I^sSz/g. To include 
this constraint in the summation, and to transform n! to a more tractable form, an identity 



n\ 



V / dt{Vtye-''\ (6) 







is used, and n in (Kt)" is replaced by Y.s^^s- Using Eqs.(4) and (6) in Eq.(5), we can 
combine K^" with {yt)"^. Summing it over Vg, we get. 



Z(^) 



/•CO 1 

V dtexp{-Vt-Y.-K'M), (7) 
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ZoZre 

where Kg{t) is a redefinition of tlie Kg by replacing 1/V by t in tlie riglit-liand side of 
Eq.(4). Summing Kg{t)/s over an integer s from 1 to oo using an identity, 

ln(l + a;) = -E^, (8) 

m=l '"' 

and carrying out the summation in Kg{t) over an even integer / including zero (Bose 
statistics) by the use of an identity, 

we obtain a final form, 

gt 1 



Zifi) = ZoZre n dtexp(-Vt) TT I 1 + 

^0 p=0 V 



f3 (ep - fiy 



/ sinh/Jy^(e,-/iP + g (,^ _ ^) 

i sinh/?(e,-^) ./(e, - ^)2 + f ' ■ ^ ^ 



In view of Eq.(lO), the interaction does not produce a new pole in Z{fi), the conventional 
definition of the BEC Instead, we obtain a new insight on the zero of Z{fi), the definition 
of the GLC. 

To extract qualitative features of Z{fi) from Eq.(lO), a most important part of the 
integrand is (l + 4 / ^ p J , which comes from / = 0, and, because of (e^ — /i)^ > /x^, 
especially important is its least term (1 + g-r) coming from a zero momentum. When 
the negative n increases, this term approaches zero first among many terms in the product 
because oi g < 0. Hence, at a critical value fic{< 0) , the zero of Z{fi) is realized by a 
cancellation of two integrals which are obtained by dividing the integrand into two parts at 
this term. Generally the density p derived by Eq.(3) is a monotonic increasing function of 
/i, to which function the zero of this Z{fi) adds a discontinuity at ;U = fie. Equations (2) 
and (3) together give the isotherm exhibiting the GLC. 
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It must be noted that the Bosons with the zero momentum play a role of a trigger to 
the GLC. The negative divergence of its p in Eq.(3) implies that the Bosons with the zero 
momentum escape from the dilute gas, thereby making a droplet . This high density assembly 
of the zero momentum Bosons is in favor of the BEC. 

Let us estimate a condition of this GLC for a weak attractive force. As a first approxi- 
mation, we can estimate pc by, 

/;.*exp(-yo(l + f^)=0. (11) 



and get a condition:/ic = —J\g\kBT /V . Further, for iJi{T)^ we use the formula in the free 
Bosegas: p{T) = -{gz/2{l)/2^fkBTBEc[{T/TBEcf'^-l]\ where (73/2(x) = Y.mX'^/nv'l^ 



Solving -^J\g\kBT/V = /i(T), we get, 

which shows non existence of the threshold strength of the attractive force: Under an 
infinite small attractive force, the GLC occurs before the chemical potential reaches zero . 

For a weak attractive force, we extend n\^ = (73/2(1), the well known condition of the 
BEC, to a general condition of the instability caused by Z{fi). {n is a number density and 
A = (27rmfcBT//i^)^^/^, a thermal wave length ||13|-) Thus we replace /x = by /ic such that 



^^^ = 5'3/2(e^^'')- Because of e^^" < 1, Tglc > Tbec for a constant n, and uglc < ^bec for 
a constant T. Figure. 2 illustrates the condition of the GLC by a solid curve for g/V = —5nK 
and by a dotted curve for g/V = —5pK with T{nK) = 30X^^ (fim) in Rb atom. (The shaded 
area represents the BEC phase of the free Bose gas: nX^ > (73/2(1).) When the temperature 
decreases and the density increases, the GLC occurs before the conventional BEC condition is 
satisfied . Under the same strength of the attractive force, it manifests itself more evidently 
in the lower temperature region. As the strength increases, its position in the phase diagram 
moves to the higher temperature and lower density region, a more realizable environment. 
Note that the GLC is very sensitive to g/V due to its cooperative nature. 

The trapped atomic gas is the Boson system in the harmonic potential. The BEC in 



the confined space has been studied theoretically from old days. Although some differences 
from that in the infinite space were reported, they do not change the essential features of 
the BEG. Hence, for the trapped Bosons as well, we can expect an essentially similar phase 



diagram [P 



The BEG with the attractive interaction was first observed in ^Li gas as a metastable 
state [|l^. A peculiar feature of the atomic gas is that we can control the sign of the 
interaction between the atoms using so-called "Feshbach resonance" . Using this technique, 
the condensate over a wide range of interaction strength is realized . By a sudden switching 
of the interaction from the repulsive to the attractive one, a collapse of the BEG gas to a 
high density spot was predicted. Recently such a phenomenon was really observed in ^^Rb 
p!6| . When such an experiment is performed in the region of nX^ < (73/2(1) (the normal 



phase of the free Bose gas) with the number of the atoms large enough to reach the thermal 
equilibrium, this type of experiment will provide us a method for finding the GLG line 
predicted in Fig. 2. 

In the attractive Bose gas, a possibility of the Boson pair has been explored by many 
people in analogy with the Gooper pair |T^. In view of the result of this paper, however, the 
attractive force responsible for such a Boson pair must have a property such that there is 
no residual attractive interaction between different Boson pairs. When such a residual force 
is there, the GLG must occur instead of the Boson pair formation. Forbidding the residual 
force, however, should impose a somewhat artificial constraint on the attractive interaction. 

This paper deals with only an initial stage of the instability. Once the Bose gas changes 
to the hquid, the short range repulsive force due to Hre plays a dominant role. (When 
the Bosons are in a high density gas state at high temperature, the situation is similar.) 
This makes Fig. 2 a more complex phase diagram: The critical point [uc, Ac) will appear on 
the GLG line, and the hne below the (ric, Ac) must be replaced by the gas-liquid coexisting 
region. Inclusion of the repulsive force will reveal a more complex nature of the later stage 
of this instability. To understand whether the formation of the droplet immediately leads 
to the BEG or not, we must get the Z{fi) which includes the Hre explicitly. The later stage 
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of this instability is an open problem. 
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FIGURES 
FIG. 1. (a) A sum of the ring diagram which represents local changes, (b) A correspond- 
ing connected diagram which represents macroscopic changes, (c) Examples of building block 
Ks{s = 1,2,3,4) of the connected diagram, which are classified by s where 2s is a number of 
particle lines. 

FIG. 2. (n, A) phase diagram of the attractive Bosons, where n is a number density and 
A = {2TTmkBT/h'^y^/^. A solid curve is a GLC line defined by nA^ = 53/2(6^^'') for g/V = -5nK, 
and a dotted curve for g/V = —5fiK with T{nK) = 30A^^(/im) in Rb atom. Shaded area is the 
BEG phase of the free Bose gas defined by nA^ > 53/2 (!)• 
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